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The lowest-lying states of the Borromean nucleus ^^Ne {^^O+p+p) and its mirror nucleus ^^N 
(^^N+n+n) are compared by using the hyperspheric adiabatic expansion. Three-body resonances 
are computed by use of the complex scaling method. The measured size of ^^O and the low-lying 
resonances of ^^F (""^^O-l-p) are first used as constraints to determine both central and spin-dependent 
two-body interactions. The interaction obtained reproduces relatively accurately both experimen- 
tal three-body spectra. The Thomas-Ehrman shifts, involving excitation energy differences, are 
computed and found to be less than 3% of the total Coulomb energy shift for all states. 

PACS numbers: 21.45.-l-v, 27.20.-|-n, 21.10.Sf, 21.10.Dr 



I. INTRODUCTION 

Nuclear halos are expected along the driplines where 
nucleon single particle s orp-states occur with sufficiently 
small separation energy [ll |31 ■ A cluster description of 
the nucleus can then be appropriate, being able to vi- 
sualize the nucleus as a system made by an ordinary 
nucleus (the core) surrounded by one or two nucleons. 
Since the Coulomb interaction works against halo forma- 
tion the appearance of proton halos requires a relatively 
light core [3, ■ Still the degrees of freedom describing 
the core and the surrounding protons may decouple and 
justify a few-body treatment. 

The lightest Borromean proton dripline nucleus is ^^Ne 
+ p) which has an odd core mass implying finite 
core spin. The low-energy properties of the two-body 
proton-core subsystems produce the two sets of known 
spin-split pairs of resonances [3- Deahng with the de- 
tails of such systems is delicate but analogous to the 
proper treatment of ^^Li |^. Unfortunately the struc- 
ture of ^^Ne, even within three-body models, seems to be 
very controversial and differing rather strongly in avail- 
able pubhcations 0,I1|- 

Nevertheless ^''Ne was recently discussed as an ex- 
ample revealing new features of the Thomas-Ehrman 
shift jlOl | . This necessarily introduces the mirror nucleus 
with less Coulomb repulsion which then must be more 
bound and possibly with a different structure. In fact 
the basic assumptions of the three-body model could be 
violated. Still it is interesting to push model applica- 
tions to test its limits. We shall therefore try to describe 
^"^N with precisely the established model parameters and 
compare effects of the (lack of) Coulomb interactions. 

The paper is organized as follows: In section II we 
describe very briefly the method used to compute three- 
body bound wave functions as well as three-body reso- 



nances. We continue in section III with the description of 
the nucleon-nucleon and nucleon-core interactions needed 
to study -'^'^Ne and "'^^N. We then in section IV compare 
the properties of the low-lying states of these two mirror 
nuclei. In section V we discuss in details the three-body 
Thomas-Ehrman shifts for all these states. We close the 
paper with a summary and conclusions. 



II. 



THREE-BODY METHOD 



The three-body wave functions are computed using the 
hyperspheric adiabatic expansion method 0, that 
solves the Faddeev equations in coordinate space. The 
wave function, 4" = X^V'i? is then written as a sum of 
three Faddeev components 'ip''^\xi,yi) (i=l,2,3), where 
{xi,yi} are the three sets of Jacobi coordinates defined 
for instance in We then introduce the hyperspheric 
coordinates, {p — \/ + y^, OLi — arctan(xi/?/i), il^,. , 
and . ) , and for each value of p we expand each compo- 
nent ^'■'^ in terms of a complete set of angular functions: 
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(1) 

When this expansion is introduced in the Faddeev 
equations they can be separated into an angular and a 
radial part. The angular part takes the form 

+ + 4P + e-^) = A„(p)<^« 

(2) 

where Vjk is the two-body interaction between particles 
j and fc, is an angular operator and m is the nor- 
malization mass. The complete set of angular functions 
used in Q are the eigenvectors of the angular part of the 
Faddeev equations, that are labeled with the index n and 
whose eigenvalues are denoted by A„(p). 

Finally, the coefficients fnip) in the expansion are 
obtained after solving the coupled set of equations given 
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by the radial parts of the Faddeev equations: 

Up) 

fn'{p)^0 (3) 
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where V^b is a three-body potential used for fine-tuning 
and the functions P„„' and Q„n' can be found for in- 
stance in TUl . The eigenvalues A„ (p) are essential in the 
diagonal part of the effective radial potentials: 



A„(p) + 15/4 
^effW = ^ 15 + ^3h(p) 



P^ 



(4) 



The hyperspheric adiabatic expansion method was ini- 
tially designed to compute three-body bound state wave 
functions, and therefore the coupled set of radial equa- 
tions (PJl was solved for radial solutions falling off expo- 
nentially at large distances. In principle the method can 
also be used to calculate continuum and resonance wave 
functions. Then we must require the correct asymptotic 
behaviour for the solutions to eqs.© as described in [T^ . 

Calculations of resonance wave functions are in prac- 
tice significantly simplified by using the complex scal- 
ing method |l3| where the radial coordinates are ro- 
tated into the complex plane by an arbitrary angle 9. 
This transformation of the Jacobi coordinates {x — > xe^^ , 
y — > ye*^) implies that only the hyperradius p is trans- 
formed {p — > pe*^), while the hyperangles remain un- 
changed. It is then known that as soon as the scaling 
angle is larger than the argument of the resonance, 
then the complex rotated resonance wave function falls 
off exponentially, exactly as a bound state. Therefore, 
after complex scaling, the resonances can be computed 
with the same numerical techniques as for bound states. 



III. 



TWO-BODY POTENTIALS 



In this section we investigate the two-body poten- 
tials needed to compute the three-body system ^^Ne 
The short-range interaction for the mirror 
nucleus ^^N (^^N+n+n) is then in principle the same 
although the assumptions of the three-body model are 
much less convincing due to the larger binding energy. 
The spin-dependence of the effective two-bod y in terac- 
tions must be carefully chosen as shown in For 
symmetry reasons the spin-spin and spin-orbit operators 
in the nucleon-nucleon interaction should be Si -82 and 
£ ■ (si -|- S2), respectively, where Si and S2 are the spins 
of the two nucleons and £ is their relative orbital angular 
momentum. 

For the nucleon-core interaction it is necessary to intro- 
duce operators that conserve the usual mean field quan- 
tum numbers, i.e. the nucleon total angular momentum 
jn = £ -\- Sn and the total two-body angular momentum 
j = jn + Sc, where £ now is the relative nucleon-core or- 
bital angular momentum, and Sc and s„ are the core and 



nucleon spin, respectively. This almost uniquely deter- 
mines the spin operators as the usual fine and hyperfine 
terms £ ■ s„ and Sc ■ jn . 



A. Nucleon-nucleon interaction 

For the nucleon-nucleon short-range interaction we use 
the operators mentioned above, and in particular the po- 
tential given in 

VNNir) = 37.05e-('-/i-3i)' - 7.38e-('^/i-«4'' 

-23.77e-(''/i-45)'.g • s + 7.16e-('^/2.43)^^^^ (5) 

(eZ„)2 



+ (49.40e"('-/i-3i)' + 29.53e"('-/i-»4)') 



Si- S2 + 



where s = Si -I- S2 and 6*12 is the usual tensor operator, e 
is the unit electric charge and Zn{— 0, 1) is the nucleon 
charge number. The strengths are given in MeV and the 
ranges in fm. This potential reproduces the experimental 
scattering lengths and effective ranges of the ^Sq, ^Pq, 
•^Pi, and ^P2 waves. We use the same interaction for 
relative orbital angular momenta larger than 1. 



B. Nucleon-core interaction 

For the nucleon-core interaction we construct an 
dependent potential of the form: 



-Erf(r/6e) , (6) 



'Sii^~ni\r)£-sn + 
r dr 



where s„ and Sc are the spin of the nucleon and the core, 
respectively, £ is the relative orbital angular momentum 
between the two particles, jn = £ + s„, and Zc is the 
proton number of the core. The error function Erf de- 
scribes the nucleon-core Coulomb interaction of a gaus- 
sian core-charge distribution where be = 2.16 fm is fitted 
to reproduce a rms charge radius in ^^O of 2.65 fm. This 
value is obtained from the measured rms charge radius 
in (2.71 fm [11) by rescahng it by an A^/^ factor. 

As discussed in [3 the choice of these spin operators 
permits a clear energy separation of the usual mean-field 
spin-orbit partners ^£+1/2 and (1-1/2- In this way it is 
possible to use a nucleon-core interaction such that the 
low-lying states have well defined ij^^ quantum numbers, 
like the pi/2 states in ^"Li or ^5/2 states in ^^F. The use 
of the Sc-Sn spin- spin and £-(sc + Sn) spin-orbit operators 
makes this impossible, since then jn is not a conserved 
quantum number and the states in the two-body system 
are necessarily mixtures oi £i^i/2 and ii-i/2 components. 
This is especially problematic in the case that one of these 
states is forbidden by the Pauli principle, like for instance 
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TABLE L Range (be) and strengths of the central (Ss*'), spin- 
spin {Si'J), and spin-orbit (Sio^) potentials in eq.Q. The 
diffuseness a is 0.65 fm in all the cases. 



£ 


be (fm) 


S^/^ (MeV) 


sis' (MeV) 


Sii' (MeV-fm^) 





3.00 


-53.91 


0.92 




1 


2.70 


-19.99 


0.69 


-25.0 


1 


2.92 


-54.15 


0.35 


-25.0 


2 


2.85 


-58.45 


0.24 


-25.0 



The shapes of the central (/c ), spin-spin (/is ) and 
spin-orbit (fsc)) radial potentials in eq.© are chosen 
to be Woods-Saxon functions, 1/(1 -I- exp((r — be)/a)), 
with the same diffuseness a in all cases. Once the range 
bi of each radial potential is chosen, the strengths Sc^\ 

Sss and Sso are adjusted to reproduce the experimen- 
tal spectrum of ^^N (-'^^N-l-n). For s-waves the strengths 

Sc^^ and si!s' are used to fit the energies of the Sy~^^ 
and the s^''^""'^' states (0~ and 1^ states). For d- waves 

(2) 

the strength Sso provides an appropriate spin-orbit split- 

(2) (2) 

ting of the ^3/2 and the dc^/2 states while Sc and Sss 
are used to reproduce the experimental binding energies 
of the d!"^j^^ and the d!"^j^^ states {2~ and 3~ states). 

The role of the spin-orbit interaction is here only to 
place the i2-siaX.es relatively high (they must remain 
unbound), and the precise energy of these states is not 
very relevant. In any case a appropriate estimation of the 
strength for the spin-orbit interaction requires knowledge 
of the d^^2 '^3/2 energies. In °N there are two 

unbound 1^/2^ doublets that could correspond to these 
states. Their experimental decay energies [1^1 (1.90 MeV 
and 2.58 MeV, or 2.27 MeV and 2.86 MeV) are used to 
estimate the strength of the spin-orbit interaction for d- 
waves. 

The value of the range parameter hi is determined by 
the fact that by switching on the Coulomb potential the 
experimental spectrum of ^^F {^^0-\-p) should be repro- 
duced. In table U we give the resulting values of the pa- 
rameters used for the Woods-Saxon radial form factors 
in eq. © . The partial waves with t = and 1 — 2 are by 
far the most important in the present context. 

The s-wave potential has a deeply bound state at —31.0 
MeV in ^^N and at -26.2 MeV in ^^p. These states cor- 
respond to the Si/2 nucleon states occupied in the ^'^N 
or the ^^O core. They are then forbidden by the Pauli 
principle, and should be excluded from the calculation. 
This is implemented as in refs. |^ IT7i | by use of the phase 
equivalent potential which has exactly the same phase 
shifts as the initial two-body interaction for all energies, 
but the Pauli forbidden bound state is removed from the 
two-body spectrum. We then use the phase equivalent 
potential of the central part of the Woods-Saxon s-wave 
potential in table Thus the s-states actively entering 
the three-body calculations are the second states of the 



TABLE II: Four lowest states in and ^"^F obtained with 
the nuclear potential specified in table |I] For ^^F we give the 
energies and widths of the two-body resonances {Er, F). The 
experimental data are from p^. Error bars are not specified 
when they are smaller than the last digit. For unbound states 
the energies are decay energies above threshold. 



r 




Exp." 


16 p 


Exp." 




-2.37 


-2.371 


(0.53,0.02) 


(0.535, 0.040 ± 0.020) 


I- 


-2.09 


-2.094 


(0.71,0.07) 


(0.728 ± 0.006, < 0.040) 


2" 


-2.49 


-2.491 


(0.96,0.01) (0.959 ±0.005, 0.040 ± 0.030) 


3- 


-2.19 


-2.193 


(1.23,0.01) 


(1.256 ±0.004, < 0.015) 



"From ref. ll6l 



Woods-Saxon potential. For the d-states no Pauli exclu- 
sion is necessary. 

The bound states in ^^N and low-lying resonances in 
^^F with J'^—Q~, 1^, 2^, 3^ are all obtained by coupling 
Si/2 and c?5/2 with the core-spin of 1/2. The calculated 
results are in table ^ compared with the experimental 
data for these states. The procedure of fitting the nuclear 
potential to reproduce simultaneously both the ^^N and 
the ^^F spectra is apparently efficient as the data is rather 
nicely reproduced. In fact this is not possible with other 
values for h^. 

In tableQlwe also specify a p-wave interaction although 
these partial waves are expected to have only insignificant 
effects. The reason is that the lowest p-shell is fully oc- 
cupied in the core and the unoccupied P3/2 orbit is above 
the c?3/2-states and even higher than the fi/2- Never- 
theless, since the calculation will include p-wave com- 
ponents at least an estimate of the parameters for the 
corresponding interaction is desirable. We do this by 
using the knowledge of the unbound 1+ and 2+ states 
in ^^N immediately above the bound 1~ state (with ex- 
perimental decay energies 0.86 MeV and 1.03 MeV [l6j . 
respectively). These resonances must arise from the cou- 
pling of a P3/2 neutron with the spin 1/2 of the core (an 
/? /2 neutron can not couple to 1 or 2) or perhaps by core 
excitation of api/2 neutron. Again the potential param- 
eters must be such that after switching on the Coulomb 
interaction the experimental decay energy of 4.30 MeV 
for the 1+ state in ^^F has to be also reproduced (the 
experimental decay energy of the 2+ state of ^^F is not 
available). The parameters fulfilling these conditions are 

given in the second line of tabled The value of S'io^ has 
been arbitrarily chosen to be the same as for d-waves. 

The lowest p shell is fully occupied in the ^^N or -'^^O 
core. We should then apply the same treatment as for 
s-waves to the p-wave nucleon-core interaction, using a 
potential with deeply bound states that are afterwards re- 
moved by the corresponding phase equivalent potentials. 
For consistency we also tested a deep p-wave potential as 
given in table ^ with range and strengths comparable to 
the s and d potentials. The bound states in these deep 
s and p potentials produce a charge distribution with a 
rms radius in -'^^O of 2.63 fm consistent with the value 
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used in the Coulomb potential. Furthermore the bind- 
ing energies of the and P3/2 states in ^'^O are —7.29 
MeV and —11.3 MeV, resp ectively, both consistent with 
the experimental data [l8| . 

Nevertheless, since the p waves basically have no effects 
in the three-body calculation we use for simplification 
the shallow i—1 potential given in table |l] without bound 
states. In this way the computing time is significantly 
reduced without loss in the computations accuracy. 



IV. RESULTS FOR "NE AND ^^N 

We use the two-body interactions determined as de- 
scribed in the previous section. The low-lying nucleon- 
core valence space is expected to consist of s and d waves. 
With spin and parity of 1/2 ~ for the core and two identi- 
cal nucleons in the sd valence space we can construct to- 
tal angular momentum and parity states with J'^=l/2~, 
3/2", 5/2-, 7/2-, and 9/2-. The next shells {Jy2, P3/2, 
• ■ • ) may also contribute but significant amounts of such 
components also indicate similar contributions from core 
excitations. These structures involve particle-hole exci- 
tations either from the sd to the p/-shell or from the p 
to the sd-shell. We shall neglect these core excitations. 



A. Components 

To solve the eigenvalue problem given in eq.Q 
we expand the angular eigenvectors in the basis 
{yttyA'^^^^^i^^Vr) ®Xs^sy,s}, where yf^i^,L are the 
hyperspheric harmonics and x is the spin function |ll| . 
For each of the three Jacobi coordinate sets i the coor- 
dinate Xi is the vector connecting particles j and fc, the 
quantum number is the relative orbital angular mo- 



p(fm) 
10 



mentum of particles j and k, £y is the relative orbital 
angular momentum of particle i and the center of mass 
of the jk two-body system. The spin Sx is the coupled 
spin of particles j and k, and Sy is the spin of particle i. 
Finally L and S are the coupling of and £y , and of 
and Sy , respectively, and they couple to the total angular 
momentum J of the system. The hypermomentum K is 
given by 2n + + £y where n is a non-negative integer 
counting the number of nodes in the Jacobi polynomials. 

The first step in the calculation is then to choose the 
components to be included in the expansion of the angu- 
lar eigenvectors. By direct but extensive computations, 
we have found that the components needed for ^^Ne are 
essentially s, p, and d-waves. Only for high angular mo- 
mentum (J=7/2 and 9/2) higher partial waves can be 
relevant. We then use the same components for ^^N. 

After solving the angular part of the Faddeev equations 
(PI we extract the angular eigenvalues A„(p), that deter- 
mine almost entirely the effective potentials entering in 
the radial equations For both ^^Ne and ^^N we also 
here maintain the same number of lowest-lying adiabatic 
potentials for use in the radial equations Q . We compute 
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FIG. 1: The four lowest angular eigenvalues A„(p) for the 
1/2', 3/2", 5/2", 7/2" and 9/2" states of ^'^Ne as function 
of p where the normalization mass m equals the nucleon mass. 



first the bound state solutions falling off exponentially at 
large distances. Then the resonance eigenfunctions are 
found in complete analogy as exponentially falling solu- 
tions to the similar equations obtained by complex rota- 
tion of the hyperradius. 



B. Spectrum of ^^Ne 

We show the results in fig^ for the four deepest ef- 
fective potentials for the 1/2", 3/2", 5/2-, 7/2", and 
9/2" states in ^^Ne. It is known that at p—0 the values 
of the A's must reproduce the hyperspherical spectrum, 
K{K+A) [HI. In our case of positive total parity in the 
valence space K must be even, i.e. K=0,2,4,- ■ ■ . 

In the figure we observe that the A function start- 
ing at zero (K—O) does not appear. This is due to the 
phase equivalent s-wave potential between proton and 
core where the deepest state of the initial pot ential is re- 
moved to account for the Pauli principle By using 
the initial deep two-body potential instead we obtain a A 
function starting at zero for p = and diverging to —00 
at large distances. This behaviour of the lowest A char- 
acterizes the existence of a bound two-body state [lH . 
This state is actually the Pauli forbidden state which 
could have been computed and then omitted from the 
basis. Instead we suppressed the Pauli forbidden state 
by using the more consistent procedure with the phase 
equivalent potential. 

For short-range potentials it is also known that at in- 
finity the values of the A's must again follow the hy- 
perspherical spectrum [TH . However, this behaviour is 
changed for eigenvalues corresponding to unbound two- 
body states as soon as long-range interactions like the 
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TABLE III: The second and third columns give the experi- 
mental and computed bound state (l/2~) and decay energies 
(3/2", 5/2", 7/2", and 9/2") in ^^Ne (in MeV). The fourth 
column are the strengths S (in MeV) of the gaussian three- 
body forces that for a range of 4.0 fm give rise to energies 
matching the experimental values. The fifth column gives the 
expectation value of the three-body force for the correspond- 
ing ^^Ne solutions. The last column is the contribution to the 
norm of the first three terms in the expansion 



J" 




Ecomp 


S 




An = l,2,3 (%) 


1/2- 


-0.94 


-0.79 


-0.6 


-0.2 


88.5, 11.1, 0.4 


3/2" 


0.34 


0.63 


-1.4 


-0.3 


90.7, 8.9, 0.2 


5/2- 


0.82 


0.91 


-0.4 


-0.1 


77.2, 16.9, 5.6 


7/2- 


2.05 


2.24 


-0.8 


-0.2 


97.5, 2.3, 0.2 


9/2- 


2.60 


2.70 


-0.1 


-0.1 


91.9, 4.2, 3.8 



Coulomb potential are present. The reason is that the 
influence of the short-range interactions then disappear 
outside p- values corresponding to a few times the range 
of the interaction whereas the Coulomb potentials mul- 
tiplied by give rise to linearly increasing A functions 
even at asymptotically large distances 0|. This linear 
increase must appear as soon as only the Coulomb poten- 
tial has an influence. The slopes depend on the geometric 
structure of the three-body system as the size increases. 

The ground state of ^^Ne is bound, and has a two- 
proton separation energy of —944 keV. The structure 
is about equal amounts of proton-core and c?^-waves. 
The computed root mean square radius is 2.8 fm consis- 
tent with the experimental value of 2.75 ± 0.07 fm [20l |. 
All the excited states are unbound and computed by ap- 
plication of the complex scaling method. The excitation 
energies of the two lowest excited states are 1288±8 keV 
for the 3/2" state and 1764±12 keV for the 5/2" state 
For both these states proton-core s-d mixed com- 
ponents are dominating. In |^ also a 7/2- and a 9/2" 
states are reported with excitation energies 2997±11 keV 
and 3548±20 keV, respectively. These four excitation en- 
ergies correspond to the decay energies (energies above 
threshold) given in the second column of table UTTl For 
these states the c?^-waves dominate. More details about 
the structure is available in j21| . 

The resonances obtained for -'^'^Ne are extremely nar- 
row with widths much smaller than the accuracy of our 
calculations. Thus, application of the complex scaling 
method allows the use of very small scaling angles. Typ- 
ically complex scaling angles of 6*^10-^ are able to find 
the ^^Ne resonances. For these scaling angles the com- 
plex scaled A's can hardly be distinguished from the non- 
rotated functions in fign The imaginary parts are very 
small and would appear on the zero line if plotted on the 
figure. 

Using the A functions in figQ] we obtain the ^^Ne 
ground state binding energy and the decay energy of the 
excited states shown in the third column of table IIIII 
As seen in the table the computed states are systemati- 
cally less bound than the experimental value. This fact 



is actually expected, since three-body calculations using 
pure two-body interactions typically undcrbind the sys- 
tem. This problem is solved by inclusion of the weak 
effective three-body potential V^b in that accounts 
for three-body polarization effects arising when the three 
particles all are close to each other. Therefore the three- 
body potential has to be of short range, while the three- 
body structure essentially is independent of the precise 
shape. This construction furthermore ensures that the 
two-body resonances remain unaffected within the three- 
body system after this necessary fine-tuning. The effec- 
tive total potential entering is then given by Q). 

The precise range of the three-body interaction also 
plays a limited role. This is because the three-body force 
is very weak compared to the depth of the full potential 
and furthermore it is largest for small p- values, where the 
total potential is highly repulsive. It is then clear that 
the main structure of the system can not be significantly 
modified by the choice of one or another of such three- 
body interactions. In table UTTl we give the strengths of 
the gaussian three-body potentials which for range equal 
to 4 fm are needed to match the experimental energies of 
all these (ground and) excited states. One way to mea- 
sure and compare the effect of the three-body force in the 
different calculations is to compute the expectation value 
(V3fc(p)) of the three-body potential for the correspond- 
ing J'" solutions. In table ITTTl we also give this quantity 
which measures the contribution of the three-body force 
to the energy of the three-body system. A variation of 
the range of the three-body force within reasonable limits 
is not modifying the results. 

In the last column of table IIIII we give the contribu- 
tion to the norm of the wave function of the first three 
terms in the expansion Typically only two terms 

are enough to get an accuracy of 99%, and only for the 
5/2" and 9/2" states the third term is giving a sizable 
contribution. 

The spectrum of ^''Ne has been previously investigated 
in 0, 111- In both works the 3/2" and 5/2" levels are 
reversed compared to the experimental data, although 
in this deficiency is corrected by use of an appropri- 
ate three-body interaction. In the present work these 
problems are not encountered. When only the two-body 
forces describing properly the ^^F spectrum are used, the 
ordering in the computed ^''Ne spectrum is correct, as 
seen in the third column of table IIIII Then, the use of 
a small effective three-body force is enough to fit the ex- 
perimental data. 



C. Spectrum of "N 

Interchanging all neutrons and protons in ^^Ne leads 
to the mirror system -'^'^N which then analogously should 
be described as a three-body system with the ^^N-core 
surrounded by two neutrons. The structure should then 
be obtained simply by switching off the Coulomb inter- 
action for ^''Ne, since the strong interaction is precisely 
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TABLE IV: ^^N spectrum obtained using the same nuclear 
two-body interactions as for ^'^Ne, and the effective three- 
body force given in table lllll The third column gives the 
results obtained with the Argonne (A) nucleon-nucleon po- 
tential plus the Woods-Saxon (WS) nucleon-core interaction 
(table 13 • The last column gives the experimental data ll^l- 



p(fm) 

10 



r 


E (MeV) 


Sa+ws (MeV) 


Exper. 


1/2- 


-8.54 


-8.31 


-8.374 




-3.72 


-3.66 


-4.711 


3/2" 


-6.63 


-6.80 


-7.000 




-3.83 


-5.03 


-5.174 


5/2" 


-6.32 


-6.36 


-6.467 


7/2" 


-5.24 


-5.17 


-5.245 


9/2" 


-4.58 


-4.59 


-4.745 



the same due to charge symmetry. In this way we can 
compute the properties of ^^N. 

The results are listed in column two of table IIVI and 
not surprisingly stronger binding is obtained. First, ^''N 
is not a Borromean system. The number of three-body 
bound states has also increased substantially, i.e. we find 
two bound states both for l/2~ and 3/2~, and one for 
5/2~, 7/2~ and 9/2~. The computed energies of these 
states agree pretty well with the experimental values |16| . 
The discrepancy with the experiment is always smaller 
than 5%. Only for the excited l/2~ and 3/2~ states a 
larger disagreement with the measured energies appears. 
In the third column of the table the calculations labeled 
by A+WS use the accurate Argonne nucleon-nucleon 
potential denoted in by vis, and the Woods-Saxon 
nucleon-core potential in tabled 

The relatively small differences between the mirror nu- 
clei are perhaps not as self-evident if we consider the be- 
havior of the angular eigenvalues shown in fig|21 We plot 
only the three lowest functions used to compute the dif- 
ferent states in ^^N. As for -'^^Ne there is no A function 
starting at zero due to the removal of the Pauli forbid- 
den s-state by use of the phase equivalent potential. The 
main difference compared to figO is the divergence to- 
wards — oo of all the A functions. This is a reflection of 
corresponding bound states in the two-body subsystems 
consistent with the quantum numbers of the three-body 
system. This underlines that ^^N cannot be a Borromean 
nucleus TTl. The qualitatively different behavior seen 
in figsQ] and 121 also emphasizes that the agreement with 
measurements for both nuclei is not a trivially build-in 
property of the present description. The model is consis- 
tent in a more profound way. 

The spatial structure of ^^N is less extended than that 
of ^^Ne because its larger binding energy. The length 
scale as defined in 2.3j is po ~ 5.5 fm (po « 5 fm for 
^^Ne) and the dimensionless measures of size and bind- 
ing energy are 

{{p/Pof) ~ 0.9 and mBpl/h^ « 6.0. The 
ground state of ^^N is located in the same region as ordi- 
nary nuclei. Among the excited states shown in table Hvl 
the less bound is the second l/2~ state. For this state the 




10 10 10 20 

P (fm) P (fm) p (fm) 



FIG. 2; The three lowest effective potentials A„(p) for the 
1/2" and 3/2" (upper part), and 5/2", 7/2", and 9/2" 
(lower part) states of ^^N. 




^"'^e (fm) 



,,o-P^' 



FIG. 3: Contour diagram for the probability distribution 
of the 1/2" ground state of ^'^N. The square of the three- 
body wave function is integrated over the directions of the 
two Jacobi coordinates. 



corresponding dimensionless size and binding are 1.9 and 
3.4, respectively, still located in the region of ordinary 
nuclei. The three-body structure is further illustrated in 
figEl where we show the square of the three-body wave 
function integrated over the directions of the Jacobi co- 
ordinates and multiplied by the volume element. The 
structure resembles that of ^^Nc with two similar domi- 
nating peaks. 

The known properties of ^^N seem to be rather well 
reproduced with the model parameters for ^^Ne. An ad- 
ditional very small retuning of the strength of the effec- 
tive three-body interaction could fit the lowest of each 
J'^ states. However, this would not improve the agree- 
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ment of the second l/2~ and 3/2~ states which then 
move into the positions —3.62 MeV and —4.14 MeV, re- 
spectively. This is probably because other effects are im- 
portant, e.g. different components could now contribute 
both from valence space and from excitations. This is 
equivalent to an attempt to describe higher-lying reso- 
nances in ^''Ne. They may also require another and per- 
haps enlarged Hilbert space. Further investigations of 
these well bound excited states are beyond the scope of 
the present report. 



V. THOMAS-EHRMAN SHIFTS 

The only difference in the computations of the two 
mirror nuclei is omission of the Coulomb interaction for 
^''N. This similarity is assumed inherently to describe the 
fundamental charge symmetry of the strong interaction. 
The immediate implication is that the differences in the 
spectra entirely must be produced by the Coulomb po- 
tential. The obvious difference is the shift of all energies 
towards stronger binding when the Coulomb potential 
is suppressed. However, this trivial overall shift is ac- 
companied by a modified structure of the states. This is 
especially seen for the s-wave components which are less 
influenced by centrifugal barrier effects. The Coulomb re- 
pulsion tends to increase the size of a given state simply 
by minimizing the energy. The s-waves are here more 
influenced than higher partial waves and the shifts are 
then larger. This double difference in energy (excitation 
energy difference) is for single-particle energies called the 
Thomas- Ehrman shift It is due to the Coulomb in- 
teraction but not necessarily in a straightforward way. In 
a recent work the Thomas-Ehrman shift was investi- 
gated in the three-body systems ^^O and -'^^Ne. 

We compare in tabled the spectrum of excitation en- 
ergies for the two mirror nuclei. The second and third 
columns give the experimental excitation energies for the 
different states in ^^Ne and ^^N, respectively. The gaus- 
sian nucleon-nucleon interaction in eq. ^ and the Woods- 
Saxon nucleon-core potential in table together with the 
three-body forces in table lllll reproduce the experimen- 
tal -'^'^Ne excitation energies. This calculation is denoted 
by G-I-WS. As in table ^3 the calculation denoted by 
A-I-WS uses the Argonne nucleon-nucleon potential and 
the Woods-Saxon nucleon-core potential in table P A 
small three-body force also permits to reproduce the ex- 
perimental ^^Ne spectrum. 

When these interactions are used for ^^N the energies 
in table IIVI and in the fourth column of table |Vl are 
obtained. The computed excitation energies (^'^'N)th of 
-'^^N are systematically higher than those measured. Still 
the agreement is surprisingly good in view of the fact 
that each state is computed independently by expansion 
on individual basis components without any parameter 
adjustment. This agreement is especially good for the 
A-t-WS calculation, where the short distance properties 
of the nucleon-nucleon interaction are carefully treated. 



TABLE V: The first column indicates the different states 
and calculations performed. The initials refer to the nucleon- 
nucleon-f-nucleon-core interactions used (G: Gaussian, WS: 
Woods-Saxon, A; Argonne). The second column gives the 
measured excitation energies of low-lying states in ^''Ne. The 
computed ^^Ne spectrum is identical. The measured and com- 
puted (^^N)th spectra in the mirror nucleus ^'^N are given in 
the third and fourth columns. The fifth and sixth columns 
contain the experimental (Ae) and computed (Ath) Thomas- 
Ehrman shifts. The seventh column shows the difference 
Ae-th be twee n experimental and computed energies for ^^N 
(see table ITVl . All the energies are given in MeV. 



J" 


i^Ne 






A, 


Ath 


Ae-th 


1/2- 


G-hWS 


0.0 


0.0 


0.0 


0.0 


0.0 


0.17 


1/2- 


A+WS 




55 


55 


55 


55 


-0.06 


3/2- 


G+WS 


1.29 


1.37 


1.91 


0.08 


0.62 


-0.37 


3/2" 


A+WS 




55 


1.51 


5) 


0.22 


-0.20 


5/2" 


G-I-WS 


1.76 


1.91 


2.22 


0.15 


0.46 


-0.15 


5/2- 


A+WS 




55 


1.95 


55 


0.19 


-0.11 


7/2- 


G+WS 


3.00 


3.13 


3.30 


0.13 


0.30 


-0.01 


7/2- 


A+WS 


;5 


55 


3.14 


55 


0.14 


-0.08 


9/2- 


G+WS 


3.55 


3.63 


3.96 


0.08 


0.41 


-0.17 


9/2- 


A+WS 


;5 


55 


3.72 


55 


0.17 


-0.16 



Furthermore the states in ^^N are well bound and the 
assumptions of independent degrees of freedom in the 
three-body cluster model cannot be very well fulfilled. 

The experimental shifts (Ae) are given in the fifth col- 
umn of table The total Coulomb shift (given by the 
energy difference between a -'^''Ne state and the corre- 
sponding ^^N state) is of around 7.4 MeV. The values of 
Ae are then remarkably small compared to the 10% of 
the total Coulomb shift for the classical example of the 
single-particle s and d-states in ^^O and "'^'^F. It is tempt- 
ing to conjecture that this is due to the stronger effect 
of the centrifugal barrier in the three-body system where 
even the s-states feel a barrier. The Coulomb repulsion 
effect is then less pronounced than for a two-body system 
where the absence of the centrifugal barrier is the basic 
explanation. 

The computed shifts {Ath) are shown in the sixth col- 
umn of the table. For the G+WS calculation they are 
clearly larger than the measured values, but also in these 
cases the computation represents an accuracy better than 
10% of the Coulomb shift. In the A+WS computation 
Ath is no more than 3% of the Coulomb shift, and shows 
a better agreement with the experiment. 

However, the {^'^N)th values given in the fourth column 
of table Ivl are obtained by comparison to the computed 
ground state 1/2^ energy for each calculation. This en- 
ergy differs from the experimental value by 170 keV and 
—60 keV for the G+WS and A+WS calculations, respec- 
tively (see last column in table . These numbers do 
not enter when comparing to the total two-nucleon sep- 
aration energy, but they do enter in the computed shifts 
(Ath) in the sixth column of tabled Therefore the un- 
certainty refiected in these computed ground state ener- 
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gies is comparable to the experimental Thomas-Ehrman 
shift Ae we are trying to reproduce. A dedicated effort 
is needed to reduce these uncertainties. 

In these comparisons the measured values include all 
many-body effects while the computations are within 
the three-body model. To estimate effects of structure 
changes we can compare the properties of these mirror 
nuclei through direct computations of Coulomb energies 
with the model wave functions. Following we con- 
sider the first-order perturbative contribution to the ^^N 
energy from the Coulomb potential, i.e., 

A(i) = (vI/(i7N)|Ko«i|*("N)) (7) 

where 4' is the three-body ^^N wave function obtained 
without Coulomb interaction between core and valence 
particles and reproducing the experimental -'^^N spec- 
trum. Then the valence neutrons are substituted by pro- 
tons in precisely the same configurations arriving at an 
artificial ^^Ne wave function. Then Vcoui is the resulting 
Coulomb interaction between the three pairs of charged 



particles. Thus Ac is the diagonal contribution to the 
Coulomb shift if the wave function remains unchanged. 
In ,2^1 the difference 

Ate = Ae - A^D (8) 

is referred to as the Thomas-Ehrman shift, where Ac is 
the experimental shift between the ^^N and ^^Ne ener- 
gies. Then Ate represents the reduction in the Coulomb 
energy in ^^Ne produced by the modified structure in the 
single particle states. 

In table IVII we give these Thomas-Ehrman shifts 
(Ate) arising from experimental (Ac) and computed 
(Ai^-*) Coulomb shifts between the mirror nuclei ^^Ne 
and ^^N. We also give A5, that is an estimate of the 
Coulomb shifts due to changes of structure included in 
the three-body model. Again we give the results for the 
G-f WS and A-I-WS. AU the Ate are less than 3% of the 
diagonal Coulomb shift. 



TABLE VI: Experimental Coulomb shift (Ac), first order contribution 
(Ai^', eq.Q), Thomas-Ehrman shift (A te) and As = Ate — Ae-th for 
the 1/2", 3/2", 5/2", 7/2", and 9/2" states. G+WS and A-^WS refer 
to the nucleon-nucleon-|-nucleon-core interactions used in the different 
calculations (G: Gaussian, WS: Woods-Saxon, A: Argonne). All the 
energies are given in MeV. 



J" ^ 


1/2- 


3/2- 


5/2- 


7/2- 


9/2- 




G-I^WS A+WS 


G+WS A+WS 


G+WS A+WS 


G+WS A+WS 


G+WS A+WS 


Ac 


7.43 


7.34 


7.29 


7.30 


7.35 


A<^) 


7.64 7.53 


7.45 7.60 


7.44 7.51 


7.58 7.56 


7.53 7.53 


Ate 


-0.21 -0.10 


-0.11 -0.26 


-0.15 -0.22 


-0.28 -0.26 


-0.18 -0.18 


As 


-0.38 -0.04 


+0.26 -0.06 


0.00 -0.11 


-0.27 -0.18 


-0.01 -0.02 



The values of Ate obtained are again highly influenced 
by the structure of the "'^^N states with the different cal- 
culations. As mentioned above, the agreement between 
computed and experimental two-neutron separation en- 
ergies in ^^N can be considered rather good (see table 
IIV|) . The experimental energy is recovered for the calcu- 
lations in table lVTl bv including in each case the appropri- 
ate three-body interaction, that as we know, keeps almost 
unchanged the three-body structure. From table HVI we 
observe that in some cases the ^^N states are up to 0.2 
MeV more bound in one of the computations compared to 
the other. These states are then more compact, and the 
Coulomb repulsion should in principle be larger. This 
is clearly seen in table IVII where the larger values for 
Ai^' for the 1/2", 3/2-, 5/2", 7/2", and 9/2- states 
are respectively the G+WS, A+WS, A+WS, G+WS, 
and A+WS calculations, that are precisely the compu- 
tations giving the more bound state for each level. In 



the 9/2 case, since the binding energy with the G+WS 
and A+WS calculations es pretty much the same, then 
also the Ac^^ value is the same in both cases. 



The main conclusion after analysis of tables Ivl and IVII 
is that the computed Thomas-Ehrman shifts are highly 
determined by the detailed structure of the ^^N states 
(for instance different three-body forces can significantly 
change the results). The small change in the structure 
from calculation to calculation is in our case important 
enough to produce large uncertainties in the computed 
Thomas-Ehrman shifts. These uncertainties are probably 
much smaller for a system, which in contrast to ^^N, 
undoubtedly can be described as a three-body system. 
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VI. SUMMARY AND CONCLUSIONS 

The Borromean nucleus ^'^Ne and its non-Borromean 
mirror ^^N arc investigated in a three-body model where 
two nucleons surround cores of ^^O and ^^N, respectively. 
We employ the well tested hyperspheric adiabatic expan- 
sion method. Then the two-body interactions must first 
be determined to reproduce the properties of the two- 
body subsystems. We carefully choose a spin-dependent 
form of the nucleon-core interaction such that the orbits 
of both the core and the valence nucleons can be treated 
consistently to lowest order in the mean-field approxima- 
tion. Then we are guarantied that the fundamental as- 
sumption in the three-body model of decoupled motion 
of core and valence nucleons is fulfilled as well as possible. 
We then proceed to determine parameters of the inter- 
actions such that the lowest four resonance energies of 
"'^^F are reproduced. For this we use the Coulomb energy 
of a gaussian charge distribution of measured root mean 
square radius. The computed rms radius of the core is in 
agreement with the measured size. 

The three-body ground state and four measured ex- 
cited states of ^^Ne are then computed. The computed 
states are systematically slightly underbound compared 
to the experimental energies, but reproducing properly 
the experimental angular momentum ordering. Agree- 
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In conclusion, the three-body model describes effi- 
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computed three-body Thomas-Ehrman shifts are then 
meaningful although relatively inaccurate. 
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